Polyelectrolyte Electrophoresis in Nanochannels: A Dissipative Particle Dynamics 

Simulation 
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We present mesoscopic DPD-simulations of polyelectrolyte electrophoresis in confined nanogeome- 
tries, for varying salt concentration and surface slip conditions. Special attention is given to the 
influence of electroosmotic flow (EOF) on the migration of the polyelectrolyte. The effective poly- 
electrolyte mobility is found to depend strongly on the boundary properties, i.e., the slip length and 
the width of the electric double layer. Analytic expressions for the electroosmotic mobility and the 
total mobility are derived which are in good agreement with the numerical results. The relevant 
quantity characterizing the effect of slippage is found to be the dimensionless quantity k 5b , where 
5b is the slip length, and an effective electrostatic screening length at the channel boundaries. 

PACS numbers: 82.35.Rs 47.57.jd 47.61.-k 82.45+z 83.50.Lh 



I. INTRODUCTION 

In recent years there is growing interest in techniques 
for manipulating single nanoparticles or macromolecules 
in micro- and nanochannel systems. The flow profiles in 
these channels and the motion of the macromolecules can 
be controlled on the nanoscale by pressure gradients and 
electric fields, and by exploiting smart channel geome- 
tries. This explains the great potential and the broad 
applicability of nanochannel devices, e.g., for analyzing 
tiny DNA or protein samples by electrophoresis^"— . 

Such systems represent a challenge for theory and com- 
puter simulation due to their high complexity. In biotech- 
nological applications, the molecules of interest are often 
charged and dissolved in buffer solutions with high salt 
concentrations. Thus electrostatic and hydrodynamic ef- 
fects compete with each other, resulting, among other, in 
a remarkable 'electrohydrodynamic screening' effect: In 
free solution electrophoresis {i.e., in a constant electric 
field), the counterion layer surrounding a charged par- 
ticle not only screens the electrostatic interactions, but 
also the dominant contribution to the hydrodynamic in- 
teractions, i.e., those that are generated by the applied 
fieldi^ii^. Therefore, simulations of electrophoresis that 
neglect the electrostatic and hydrodynamic interactions 
altogether often give results that are in good semiquan- 
titative agreement with experiments^ d^d^ . More sophis- 
ticated approaches that still allow to avoid the explicit 
representation of charges have been devised as wel l^^d"^ . 
Nevertheless, it is clear that such simplified treatments 
disregard important physics, especially in confined ge- 
ometries where electrostatic interactions compete with 
the regular steric interactions with the confining walls-®. 
Simulations that take full account of electrostatic and 
hydrodynamic interactions are clearly desirable. Such 
simulations have recently been carried out for polyelec- 
trolyte electrophoresis in free solutions"'^^'^^, but, to the 
best knowledge of the present authors, not yet for mi- 
crochannels. 



In many microchannels, an additional effect comes 
into play, which significantly modifies the effective elec- 
trophoretic response of particles to electric fields: The 
same electric field that drives the polyelectrolyte may 
also induce a total net fluid flow in the microchannel, the 
'electroosmotic flow' (EOF). Many materials commonly 
used in microtechnology like PDMS (Polydimethylsilox- 
ane) acquire charges if brought in contact with water, 
either by the ionization or dissociation of surface groups 
or the adsorption of ions from solution. To screen the 
charges on the chaimel walls, a diffuse layer of oppositely 
charged ions forms in front of the walls. In the presence 
of an external electric field E^^^*\ these ions are pulled 
along, dragging the surrounding fluid with them. One 
gets a characteristic 'plug' flow profile which saturates at 
a fluid velocity 



(1) 



outside of the diffuse layer, with the so-called electroos- 
motic mobility /i^oF- The effective migration speed vp 
of particles in microchannels in response to the exter- 
nal fields then results from two contributionsi^: The 
bare electrophoretic mobility fi^, of the polyelectrolyte 
in a fiuid at rest, and the background EOF velocity, 
Vp = VjjjQp -|-/ieE''^^*). Experimentally, it has been found 
that the former may even dominate over the latter, such 
that the polyelectrolyte effectively migrates in the direc- 
tion opposite to the applied field^. 

Now if the diffuse layer is thin compared to the chan- 
nel dimensions, Eq. ([1]) can be regarded as an effective 
boundary condition for a steady-state EOF velocity field 
Vj5Qp(r) inside the channel. (We note that in a steady- 
state situation, the external field in the vicinity to a wall 
is necessary parallel to the wall.) Cummings et al.-- have 
shown that for laminar incompressible flow, this bound- 
ary condition effectively defines the flow inside the chan- 
nel. Provided Eq. ([T]) also holds at the inlet and outlet 
boundaries of the channel, it becomes valid everywhere 
in the channel^^. Due to this remarkable similitude be- 
tween the steady-state velocity field of an EOF and the 
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externally applied electric field, the net electrophoretic 
velocity of nanoparticles or polyelectrolytes can be writ- 
ten approximately as 

VP = (AiBOP+Me)E(-*) (2) 

with the effective total mobility fit — (J-eof +/^e- We note 
that this expression relies on two assumptions which are 
both not obvious: The nanoparticles or polyelectrolytes 
stay well outside the diffuse layer covering the wall, and 
they themselves do not influence the EOF. 

The EOF amplitude at planar walls with no-slip 
boundary conditions has been calculated a long time ago 
by Smoluchowski^'^. On the nanoscale, however, the no- 
slip boundary condition does not necessarily apply. Ex- 
periments have indicated^i^, that the velocity profile 
is not strictly continuous at walls, i.e., fluids exhibit a 
certain amoung of slippage. This effect can be enhanced 
significantly by using superhydrophobic walls which are 
covered by a thin gas layer. It can also be tuned to some 
extent by designing nanopatterned surfaces with alter- 
nating hydrophobic and hydrophilic sections^''. In all of 
these cases, the appropriate mesoscopic boundary condi- 
tion is the 'partial slip' boundary condition, 

d 

Vx{±zb) ^ TSB-g^Vx{z)\,^±,^, (3) 

where zb is the position of the hydrodynamic boundary 
(which is usually close to the physical boundary, but not 
necessarily identical), and the 'slip length' 5b character- 
izes the amount of slippage. No-slip boundaries corre- 
spond to Sb — 0, full-slip to Sb oo. From Eq. ([3]), 
one would expect that the EOF amplitude is enhanced 
in the presence of slippage^^, and this is indeed found in 
simulations^^i^. The effect of slippage on EOF has been 
calculated within the linearized Poisson-Boltzmann the- 
ory, the Debye-Hiickel approximation, by Joly et al.— . 
Below, we will derive a general expression which is also 
valid beyond the Poisson-Boltzmann theory. 

In this paper we present Dissipativc Particle Dynam- 
ics (DPD)-simulations of polyelectrolyte electrophoresis 
in microchannels with varying slip lengths, at varying 
salt concentrations. We treat the solvent and all ions ex- 
plicitly, and all charges interact via unscreened Coulomb 
interactions. This allows to investigate the interplay of 
solvent and polyelectrolyte, of electrostatic and hydro- 
dynamic interactions, of electrophoresis and EOF in full 
detail, with almost no approximation. (Our only approx- 
imation is to neglect image charge effects, i.e., the di- 
electric constant is taken to constant everywhere.) Our 
results indicate that the hydrodynamic boundary condi- 
tions strongly influence the total mobility of the polyelec- 
trolyte, and the total mobility can be tuned from positive 
to negative by varying the slip length. The simulation 
data are compared with a simple analytical expression, 
which is derived based on the assumption that the flow 
profile in the channel follows the Stokes equation. The 
numerical results are in very good agreement with the 
theory. 



The paper is organized as follows. The theory is pre- 
sented in section 2. Section 3 focuses on the simulation 
method and the parameters used in our simulations. The 
numerical results will be shown in section 4. We conclude 
with a brief summary in section 5. 

II. THEORETICAL CONSIDERATIONS: EOF 
IN THE PRESENCE OF SLIPPAGE 

We consider for simplicity a planar slit channel with 
identical walls at z = ±L/2, exposed to an external elec- 
tric field Ex in the x direction. The electrostatic poten- 
tial ^{x,y,z) then takes the general form ^{x,y,z) = 
i^iz) + Ex X + const, where we can set ■0(0) = for sim- 
plicity. The electrolyte in the channel is taken to contain 
n different ion species i with local number density pi{z) 
and valency Zi, which results in a net charge density 
p{z) = Sr=i(^«^)^'*(^)- '^^'^ electric field then generates 
a force density fxiz) = p{z)Ex in the fluid. Comparing 
the Poisson equation for the electrostatic potential -0, 

d^tl^jz) ^ _p{z) 

(where is the dielectric constant), with the Stokes 
equation's 

Q'^y^ (^) 

Vs = "/^(^) = -p{z)Ex (5) 

(with the shear viscosity yy^), one finds immediately 
dzzVx{z) — dzz'4'iz) (cr Ex/rjs). For symmetry reasons, 
the profiles Vx and ip must satisfy the boundary condi- 
tion dzVx\z=o — dzil'\z=o = at the center of the channel. 
This gives the relation 

Vx{z) = ^^-^V(2) + Weof> (6) 

where we have used ■0(0) = and identified the fluid 
velocity at the center of the channel with the EOF veloc- 
ity, Vx{Q) = Weof- We further define V's := iP^^zb) (for 
no-slip boundaries, -0s is the so-called Zeta-PotentiaP'^) . 
Inserting the partial-slip boundary condition for the fiow, 
Eq. ([3]), we finally obtain the following simple expression 
for the electroosmotic mobility, 

Meof = V^of/Ex ^ Meof (1 + ^b), (7) 

where we have defined the inverse 'surface screening 
length' 

and rfi is the well-known Smoluchowski result^ for the 

^ EOF 

electroosmotic mobility at sticky walls, 

pIof = -'^ripB/r]s. (9) 
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The remaining task is to determine the screening pa- 
rameter K. If the surface charges are very small and the 
ions in the liquid are uncorrelated, it can be calculated 
analytically within the linearized Debye-Hiickel theory^!. 
The Debye-Hiickel equation for the evolution of the po- 
tential "0 in an electrolyte solution reads dzzi^ = '^ijV' 
with the inverse Debye-Hiickel screening length 



(10) 



where pi^Q is the density of ions i far from the surface. It 
is solved by an exponentially decaying function, 



III. SIMULATION METHOD 
A. Dissipative Particle Dynamics 

Newtonian fluids are effectively modeled by the Dis- 
sipative Particle Dynamics (DPD) metho d^^'^^ . DPD is 
a coarse-grained, momentum-conserving simulation tech- 
nique which creates a well-defined canonical ensemble. 

The basic DPD equations are defined in terms of the 
forces on one particle, which involve two-particle interac- 
tions 



pDPD ^ V 



(14) 



2). 



(11) 



with a conservative force Fp^ 

'■J 



Inserting that in Eq. ([81), one finds k = kjj, i.e., the sur- 
face screening length is identical with the Debye screening 
length. Eq. ^ with k — kd basically corresponds to the 
result of Joly et ali^. 

Unfortunately, the range of validity of the Debye- 
Hiickel theory is limited, it breaks down already for mod- 
erate surface potentials tpB and/or for highly concen- 
trated ion solutions. Nevertheless, the exponential be- 
havior often persists even in systems where the Debye- 
Hiickel approximation is not valid. For high ion con- 
centrations, detailed studies based on integral equa- 
tions have lead to the conclusion that the Debye-Hiickel 
approximation can still be used in a wide parameter 
range, if kd is replaced by a modified effective screening 
length^"—. For high surface charges, analytical solu- 
tions are again available in the so-called 'strong coupling 
limit', where the profiles are predicted to decay expo- 
nentially with the Guy-Chapman length2S. This limit 
is very special and rarely encountered. At intermediate 
coupling regimes, the decay length must be obtained em- 
pirically, e.g., by fitting the charge distribution p{z) to an 
exponential behavior, which is characterized by the same 
exponential behavior than 'il'{z) by virtue of the Poisson 
equation. 



E 



{Z^e) pi{z) oc 



dz^ 



(12) 



Putting everything together, the total net electrophoretic 
mobility fif (Eq. of a polyelectrolyte in the channel can 
be expressed in terms of the electroosmotic mobility /ijjoF 
as 



= 1 



fJ-e 



Me 



<,k(1 + «^b)' 



(13) 



where the ratio /ieZ/i^^Qp depends only weakly on the ionic 
strength of the electrolyte and the slip length of the sur- 
face. The main effect of slippage is incorporated in the 
factor (1 -I- K 5b)~^- 



Fij — ijUijijij). 



(15) 



(where rij denotes the distance between the centers of 
particles i and j), a dissipative force F^^ 

F,^ ^ -jDPD^D{rij){fij ■ Vij)rij (16) 

with the friction coefficient 7dpd, and a random force F^ 

F"j = \/2-foi-DkBT ujR{rij)C,ijfij. (17) 

Here C,ij — Qi is a symmetric random number with zero 
mean and unit variance, and the weighting functions of 
the dissipative and the stochastic force are related by a 
fiuctuation-dissipation theorem. 



(18) 



which ensures that an equilibrium simulation samples a 
canonical ensembl o'^^'^^ . Otherwise, the weight function 
Wdpd is arbitrary and will be chosen linear here as often 
in the literature. 



^DPD V' IJ 



1 



(19) 



where denotes the cut-off radius. 



B. Tunable slip boundaries 

The hydrodynamic boundary condition at the walls is 
realized with a recently developed method^^ that allows 
to implement arbitrary partial-slip boundary conditions: 
Wc introduce an additional coordinate-dependent viscous 
force that mimicks the wall/fluid friction 



with a dissipative contribution 

Ff = -7l wl(z) (v; - ^fwau) 



(20) 



(21) 
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coupling to the relative velocity (vi—v^aii) of the particle 
with respect to the wall, and a stochastic force 

= v/27L kBTiOLiz) X^a (22) 

which satisfies the fluctuation-dissipation relation and 
thus ensures that the local equilibrium distribution is 
again a Boltzmann distribution. Here a is a = x,y, z 
and is a Gaussian distributed random variable with 
mean zero and variance one: {xi,a) = 0, {xi,aXj,p) — 
SijSajS- The viscous coupling between fluid and wall is 
achieved by the locally varying viscosity "/^^^{z) with 
ujl{z) = 1 — z/zc up to a cut-off distance Zc. The pref- 
actor 7l can be used to tune the strength of the friction 
force and hence the value of the slip length. Within this 
approach it is possible to tune the slip length Sb system- 
atically from full-slip to no-slip, and to derive an analytic 
expression for the slip length as a function of the model 
parameters'^. 

C. Simulation details 

We have studied the electrophoresis of charged poly- 
mers of length iV = 20 in electrolyte solutions, confined 
by a planar slit channel with charged walls. All particles, 
polymer, solvent and ions, are modeled explicitly. We use 
a simulation box of size {12a x 12a x lOcr) which is peri- 
odic in X- and y-direction and confined by impermeable 
walls in the z-direction. The walls repel the particles 
via a soft repulsive WCA potential^S, of range a and am- 
plitude e. (Hence the accessible channel width for the 
particles is actually — 8a). Ions and monomers repel 
each other with the same WCA potential. In addition, 
chain monomers are connected by harmonic springs 

armonic — 

7rfc(r„- - ror (23) 

with the spring constant k = 25e/a'^ and tq = l.Ocr. Neu- 
tral solvent particles have no conservative interactions 
except with the walls. 

The wall contains immobilized, negatively charged par- 
ticles at random positions. Every second monomer on the 
polyelectrolyte carries a negative charge. The solvent 
contains the positive counterions for the walls and the 
polyelectrolyte, and additional (positive and negative) 
salt ions. All charges are monovalent, and the system as a 
whole is electroneutral. In addition to their other interac- 
tions, charged particles interact via a Coulomb potential 
with the Bjerrum length A_b = / ATrerksT = l.Ocr, and 
they are exposed to an external field — — l.Oe/ecr. 
Specifically, we have studied systems with a surface 
charge density of a a — — 0.208eCT~^. The total counte- 
rion density was p — 0.06ct~'^ and the salt density varied 
between ps =0.05625, 0.0375, 0.03, 0.025, and O.OlScr-^. 

We use DPD simulations with a friction coefficient 
7dpd = 5.0cr~^(me)^/^. The density of the solvent par- 
ticles was p = 3.75cr~^, and the temperature of the sys- 
tem was T = l.Oe/ks- For these parameters, the shear 



Table I: Slip lengths 5b for different layer friction coefficients 
jL, compared with theoretical value 3% according to Ref. 38. 



7i[a-i(me)i/2] 


3B[a] 


±3b[ct] 


S%[a] 


0.1 


14.977 


1.879 


14.000 


0.25 


5.664 


0.783 


5.458 


0.5 


2.626 


0.521 


2.613 


0.75 


1.765 


0.409 


1.664 


1.0 


1.292 


0.423 


1.190 


6.1 


0.000 


0.197 


0.000 



viscosity of the DPD fiuid - as determined by fitting 
the amplitude of Plane Poiseuille fiows^ - is given by 
T]s = (1.334 ±0.003)cr-2(me)i/2. The DPD timestep was 
St = 0.01cr(m/e)i/2. 

Tunable-slip boundary conditions were used with fric- 
tion coefficients =0.1, 0.25, 0.5, 0.75, 1.0, and 
6.1a~^{mey/^ . The range of the viscous layer was 
Zc — 2.0a. Only the solvent particles interact with 
the tunable-slip boundaries. By performing Plane 
Poiseuille and Plane Couette flow simulations with the 
above given parameters, the slip length 5b and the hy- 
drodynamic boundary positions zb can be determined 
independently^. The hydrodynamic boundary position 
is found at \zb\ = (3.866 ± 0.266)cr in all simulations. 
The corresponding slip lengths are presented in Table |T] 
together with the theoretical values predicted by the an- 
alytic expression in Ref. 38. The comparison shows that 
the simulated results are in good agreement with the the- 
ory. 

The electrostatics were calculated by P3M-*0 and 
the ELC (electrostatic layer correction)-algorithm^ for 
2D + h slabwise geometries. All simulations have been 
carried out with the freely available software package 
ESPResSoi^ii^ 



IV. NUMERICAL RESULTS 

Fig. [T] shows the average ionic distributions of an- 
ions Pa and cations pc for the salt concentration ps = 
0.05625(7"'^. Here the 'cations' include the positively 
charged salt ions and the counterions of the wall and 
the polyelectrolyte, and the 'anions' only the negatively 
charged salt ions. Due to the presence of the polyelec- 
trolyte in the middle of the channel, the average cation 
density there is slightly increased. To determine the in- 
verse effective screening length k, we have thus fitted the 
following function 

Ap = Apo(e-"" + e«")+c (24) 

to the ionic difference Ap = pc — Pa- The exponential 
fit describes the data very well (black solid line in Fig.[T]). 
The fit parameters for k are listed in Table HIl along with 
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Figure 1: Distribution of the ionic difference Ap = Pc — pa 
between cations and anions in the solution (not counting 
the polyelectrolyte) for an exemplary salt concentration of 
pa ~ 0. 05625(7 ~^ and the surface charge density a a = 
— 0.208e(j~^. The black line corresponds to an exponential 
fit (Eq. (|24|l ') with an effective inverse screening length of 
K = 2.305 ± 0.025(7^^. Inset: Distribution of cations (cir- 
cles) and anions (triangles) for the same system. 



the values for the Debye-Hiickel screening parameter no 
(Eq. ([TU])). The decay lengths are overall very different 
from those predicted by the Debye-Hiickel theory. We 
conclude that the system is outside the validity region of 
the linearized Poisson-Boltzmann approximation. This 
is perhaps not surprising, given that the individual ion 
profiles (inset of Fig. [1]) at the walls deviate strongly from 
their bulk value, i.e., these deviations can hardly be con- 
sidered as small perturbations. The surface charge is too 
high. On the other hand, the electrostatic coupling con- 
stant S — 2nZ^\%(TA ^ 0.2 (Z = 1 is the valency of 
the cations) is still much smaller than unity, hence we 
are still in a 'weak coupling' regime. This is also evi- 
dent from the fact that the effective screening param- 
eter K differs strongly from the Guy Chapman length, 
/i^^ = 2'iTXBZaA = 1.31(7. 

The EOF profiles for the same salt concentration {ps — 
0.05625(7^'^) are shown in Fig. [51 The different curves 
correspond to different hydrodynamic boundary condi- 
tions (slip lengths). As expected, the flow velocity in- 
creases drastically for larger slip lengths. All curves are 



Table II: Fitted inverse screening lengths k and Debye-Hiickel 
screening parameter k_d for different salt concentrations ps 
and the fixed counterion density of p = 0.06ct~^. 



Pslo-""] 






KB [a ^] 


0.015 


1.996 


0.041 


0.98 


0.0225 


2.011 


0.049 


1.02 


0.03 


1.983 


0.041 


1.07 


0.0375 


2.182 


0.047 


1.11 


0.05625 


2.305 


0.025 


1.21 



Figure 2: Exemplary fiow profiles for a salt concentration 
Ps = 0.05625o--^ for varying slip lengths (from bottom to top: 
5b = (14.98, 5.66, 2.63, 1.77, 1.29, 0.00)cr.) The black lines are 
the theoretical predictions obtained by integrating the Stokes 
equation (Eq. ((S])) with a fitted inverse screening length of 
K = 2.305cr~\ 




Figure 3: Ratio Meof/Meof plotted against Sbk, for the differ- 
ent salt concentrations and screening lengths given in Table |T] 
and[TTl The black line is the theoretical prediction of Eq. ^ 
with slope 1 + K(5_B. Inset: Surface potential as obtained from 
p^Q-p using Eq. ^ (circles) and indepently by a test charge 
method (triangles) as a function of the salt concentration ps. 



in good agreement with the theoretical predictions, which 
were obtained by integrating the Stokes equation 
numerically with the correct partial-slip boundary con- 
ditions. In agreement with our earlier studies at zero salt 
concentration^^, we thus find that a description based on 
the Stokes equation - a continuum equation - remains 
valid even for very narrow channels. 

The flow velocity in the middle of the channel gives 
the EOF mobility. Fig. [3] compares our numerical re- 
sults for all salt concentrations and slip lengths with the 
theoretical prediction of Eq. ([7|), where p^^^^ has been 
determined by a linear regression for each salt concen- 
tration independently. We find good agreement between 
simulation data and theory. This confirms the validity of 
our theoretical result, Eq. ([T]). It also demonstrates that 
the polyelectrolyte, which was present in all simulations, 
does not perturb the EOF even in very narrow channels. 



6 



0.14 

0.12 

S 0.1 
nj 

0.08 
qI 

1 0.06 
o 

S 0.04 
0.02 














/ ~ exp(-0.5 (z/a)-) 








a = 1.51373+/- 0.03153 














4 -3 


-2-1 1 
z[o] 


2 


3 



Figure 4: Normalized monomer distribution inside the chan- 
nel for the salt concentration ps — 0.05625(t~'^. 



The values of the EOF mobility for zero slip length, 
/^EOF' '^^^ used to determine the surface potential ipB 
for the different salt concentrations via Eq. As a 

consistency check, we have also determined ipB indepen- 
dently by inserting a test charge into the ion layer at 
z — zb- The results are shown in the inset of Fig. [31 Both 
methods give identical results. The surface potential is 
found to be largely independent of the salt concentration. 

After investigating the EOF of the solvent, we dis- 
cuss the properties of the polyelectrolyte. The proba- 
bility distribution for finding a monomer at a given po- 
sition z is shown in Fig. |4] for the salt concentration 
Ps = 0.05625(7"'^. It is approximately Gaussian with 
a peak in the middle of the channel and a variance 
Var ^ 2.28cr. Thus the polyelectrolyte mainly 'senses' 
the EOF in the middle of the channel, and the details 
of the flow profiles close to the channel walls have very 
little infiuence on its net mobility: The assumption that 
the total mobility is governed by a single EOF velocity 
■^EOF (Eq- ©) is probably legitimate. 

The infiuence of the ion profiles on the chain structure 
of the polyelectrolyte can be investigated by consider- 
ing static properties like the radius of gyration ~ 
(l/2iV2)^w^^^ _^ _ > ^j^g end-to-end ra- 
dius Rl =< (-Rat - ^i)^ >'^^- The results for these 
parameters are shown in Table IIIII Both characteristic 



Table III: Radius of gyration Rg and end to end radius R^ 
for a polyelectrolyte with A'' = 20 monomers for different salt 
concentrations ps. 





RM 


RA'A 


0.015 


3.2218 ± 0.047 


10.6480 ± 0.0314 


0.0225 


3.1661 ± 0.0041 


10.2736 ± 0.0266 


0.03 


3.1486 ± 0.0451 


10.1777 ± 0.0292 


0.0375 


3.1279 ± 0.0045 


10.0819 ±0.0287 


0.05625 


3.0825 ± 0.0045 


9.8331 ± 0.0280 
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Figure 5: Ratio /-tt/j/iEOFl plotted against &b k for all slip 
lengths (Table H]) and salt concentrations (Table [Tl|. The 
black line is the theoretical prediction of Eqn. (|13ll with ab- 
solut values of l/iEorl- the limit «; — oo, the to- 
tal mobility of the polyelectrolyte is equal to the electroos- 
motic mobility /^eof- The ratio /ie/zigoF '^^^ been fitted 
to —3.778 ± 0.128. Negative values of /ii/|MEOFl indicate 
negative total mobilities of the polyelectrolyte. Inset: To- 
tal displacement of the polyelectrolytes center of mass for 
different boundary conditions and a salt concentration of 
Ps ~ 0.05625cr~^. The total mobility becomes negative for 
I Me I IMeofI- The lines correspond from top to bottom to 
the slip lengths Sb ^ (0.00, 1.292, 1.765, 2.626, 5.664, 14.98)cr. 
Thus larger slip lengths indirectly enhance the total mobility 
of the polyelectrolyte. 



lengths decrease with increasing salt concentration due to 
a more effective screening of electrostatic interactions, in 
accordance with standard theoriea^. The ratio between 
the end-to-end radii and the gyration radii is unusually 
large, which is most likely a squeezing effect due to the 
presence of the channel walls46'4''. Specific fiow-induced 
effects such as shear-induced elongation^® are probably 
less important, since the flow profile is basically constant 
inside the channel. It should be noted that the pure 
electrophoretic mobility /ig of the chain is presumably 
modified by the confinement as an indirect effect of the 
elongation. This effect has not been investigated in the 
present study. 

The total mobility of the polyelectrolyte for varying 
boundary conditions and salt concentrations is finally 
presented in Fig. [5l The theoretical prediction of Eq. 
agrees well with the numerical results with the single fit 
parameter /ie//i°Qp — —3.778 ± 0.128. It is remarkable 
that this parameter can be set to a constant, i.e., it seems 
to be largely independent of the salt concentration ps. 
Since p^^^ does not depend on ps (see Fig. [31 inset), this 
means that pe is also independent of ps for the range of 
salt concentrations considered here^^. 

For no-slip boundary conditions with 6b ~ 0, we 
find ordinary behaviour where the polyelectrolyte fol- 
lows the applied electric field. In the presence of wall 
slip, however, the EOF becomes stronger and eventu- 
ally dominates. Then the total mobility may become 
negative, i.e., the polyelectrolyte migrates in a direc- 
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tion which is opposite to the apphed force. The inset of 
Fig. [5] illustrates this by showing the total displacement 
of the chain's center of mass for the salt concentration 
Ps = 0.05625(7"^ and various slip lengths. In nearly all 
cases except (5b ~ 0, the total mobility of the polyelec- 
trolyte is negative. 

To summarize this section, both the assumptions and 
the predictions of section |TT] are supported by our numer- 
ical results. The total mobility of the polyelectrolyte can 
therefore be adequately described by Eqs. ([7]) and p^ . 

V. CONCLUSIONS 

We have presented mesoscopic DPD simulations of 
polyelectrolyte electrophoresis in narrow microchannels, 
taking full account of hydrodynamic and electrostatic in- 
teractions. A particular focus was put on studying the 
effects of the hydrodynamic boundary conditions at the 
channel walls on the electroosmotic flow and on the net 
electrophorctic mobility of the polyelectrolyte. We have 
shown that they can be incorporated into a single dimen- 
sionless parameter k5b), where 5b is the slip length 
and K the (local) inverse screening length of the charge 
distribution at the wall. This was derived analytically 
and supported by our numerical data. It remained valid 
even for very narrow channels, where the chain confor- 
mations were affected by the confinement. 

We have shown that wall slip massively enhances the 
EOF and hence influences the total mobility of the poly- 
electrolyte. If the EOF mobility and the free drain- 



ing mobility oppose each other, i.e., if the effective 
charges on the polyelectrolyte and the walls have the 
same sign, the mobility may even become negative. As 
mentioned in the introduction, this effect has also been 
observed experimentally^. In the other case, where the 
sign of the charges on the polyelectrolyte and the wall are 
opposite, the main effect of slip is to enhance the total 
mobility of the polyelectrolyte. 

In summary, the total mobihty of polyelectrolytes in 
microchannels results from an interplay of electroos- 
motic, electrophorctic, electrostatic and slippage effects. 
The latter have a particularly strong influence and can 
be used to design channels with improved properties. For 
example, the characteristics of the channel walls could 
be designed to tune effective slip lengths^^ and hence 
flow velocities, which offers the possibility to optimize 
the time which is needed for polymer migration or sep- 
aration techniques. This could be an important aspect 
for future applications in microchannels or micropumps 
to accelerate measuring times. 
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